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ABSTRACT 

A theory is developed which relates the observed 6-parameter of a Lya ab- 
sorption line to the velocity-curvature of the corresponding peak in the optical 
depth fluctuation. Its relation to the traditional interpretation of b as the thermal 
broadening width is discussed. It is demonstrated that, independent of the details 
of the cosmological model, the differential b distribution has a high b asymptote of 
dN/db oc where m > 5, when we make the reasonable assumption that low- 
curvature fluctuations are statistically favored over high-curvature ones. There 
in general always exist lines much broader than the thermal width. We develop 
a linear perturbative analysis of the optical depth fluctuation, which yields a 
single-parameter prediction for the full b distribution. In addition to exhibiting 
the high velocity tail, it qualitatively explains the observed sharp low b cut-off - 
a simple reflection of the fact that high- curvature fluctuations are relatively rare. 
While the existence of the high b asymptote, which is independent of the validity 
of the linear expansion, is consistent with the observed b distribution, a detailed 
comparison of the linear prediction with six observational datasets indicates that 
higher order corrections are not negligible. The perturbative analysis nonetheless 
offers valuable insights into the dependence of the b distribution on cosmological 
parameters such as VL and the power spectrum. A key parameter is the effective 
smoothing scale of the optical depth fluctuation, which is in turn determined 
by three scales: the thermal broadening width, the baryon smoothing scale (ap- 
proximately the Jeans scale) and the observation/simulation resolution. The first 
two are determined by reionization history, but are comparable in general, while 
the third varies by about an order of magnitude in current hydrodynamic sim- 
ulations. Studies with non-resolution-dominated b distributions can be used to 
probe the reionization history of the universe. 
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1. Introduction 



It was recently pointed out by Rutledge ( fL997| ) that structures which arise naturally in 
hierarchical clustering models imply a high-velocity tail for the 6-parameter distribution of 
the Lya forest. The analysis was done using idealized filaments and pancakes. We show in 
this paper that the same qualitative (but different quantitative) behavior can be understood 
in the broader context of the statistics of peaks in the optical depth fluctuation. 

Both the high b tail as well as the sharp low b cut-off in the observed 6-parameter dis- 
tribution have been noted in the literature (e.g. Press fc Kybicki 199^ ; for recent results, see 
Hu et al. 1995| , [Lu et al. 1996| and |Kirkman fc I'ytler 1997[ ). Although there are subtle effects 
on the detection of narrow and broad lines due to finite S/N and continuum fitting ( [Rauch| 
et al. 19931 ), the ab ove features seem to be robust (see the above ref.). Recent hydrodynamic 
simulations and semi-analytical calculations reproduce the same features ( |Cen et al. 1994| , 
Zhang et al. 1995| , [Hernquist et al. 1996| , [Miralda-Escude et al. 199^ , [Zhang et al. 1996| , [Bi 
|fc Davidsen 199'^ and [Dave et al. 1997| ; but see also Paehnelt fc Steinmetz 1997] ) , although 
they do differ in details, a point to which we will return. 

In §1^, we propose relating the measured 6-parameter of a Lya absorption line to the cur- 
vature around a peak in the optical depth fluctuation. Under very general conditions, a high 
b tail of dN/db oc b~^ is predicted for the differential b distribution. The viewpoint adopted 
here is that the low column density Lya forest (A^m ~ 10^^ cm~^) is part of a fluctuating in- 



tergalactic medium, as predicted by structure formation models ([Bi et al. 1992 ; Reisenegger 
fc Miralda-Escude 1995[) . In we discuss the physical meaning of the measured b in the 
context of such models and its relationship with the traditional interpretation of b as the 
thermal (plus turbulence) broadening width. We illustrate these ideas by giving a concrete 
example in §|. A linear perturbative expansion of the optical depth fluctuation is developed, 
which yields a simple prediction for the full b distribution: in addition to showing the high 
b tail, it also implies a sharp low b cut-off similar to that found in the observed distribution. 
The single-parameter linear prediction, as opposed to the commonly used three-parameter 
truncated Gaussian, is compared with the b distributions of six datasets. Four of them can 
be satisfactorily described by the model while the other two cannot. We discuss the pos- 
sible causes, and argue the main reason is that higher order corrections are non- negligible. 



3 



Nonetheless, the hnear analysis provides us useful intuition on how the velocity structure 
of absorption lines reveals or depends on the cosmology, thermal history and resolution of 
observations/simulations. Finally we conclude in 



2. The High b Tail 

Let us expand the optical depth r as a function of velocity u around an absorption line: 

■1, 



t{u) = exp [Inr(u)] = r(Mmax) exp 



where Wmax is the velocity coordinate of the line center, and the prime denotes differentiation 
with respect to u and the second derivative is evaluated at Umax- The first derivative vanishes 
because r is at a local extremum. The reader will recognize the above expansion as none 
other than the Voigt profile: oc exp[— (m — 'Umax)^/^^]- The simple thermal profile suffices for 
our purpose, since we will focus on the low column density Lya forest. A similar expansion, 
applied to the density field in comoving space (as opposed to the optical depth field in 
velocity space here), has been used to study the column density distribution (|Gnedin 
|Hui 19961 ; |Hui et al. 19961) . 



In other words, b = ^—2/[\n r]". No assumption has been made about whether the 
given absorption line is thermally broadened or not. We will refer to it as the 6-parameter, 
or more simply b, instead of the commonly used Doppler parameter. As we will show in §^ 
the measured b agrees with the thermal b (oc VT) only if the neutral hydrogen density is 
sharply peaked in velocity space. Let us denote Inr by tl- The number density of peaks in 
Tl (or minima in transmitted flux) with a given b or t'I is: 

where iV„ is the number of peaks in tl per unit u, and P(r^, r'Ddr'i^dr'l is the probability 
at any given point along the spectrum that and t'I take the values in the prescribed 
ranges. The equality relating dNu/dr'l and P above follows from this argument: the number 
density of peaks for a particular realization of the spectrum is a sum of Dirac delta functions 
J2i^D{u ~ Ui) where i denotes all the places where is at a local maximum with a given 
t'I; close to a local maximum, t'j^{u) {u — Ui)T'l{ui); changing the independent variable of 
the delta function from u to r| and performing an ensemble average yields the above result 
(see [Bardeen et al. 19861) . 



It can be seen that as long as -P(t^ = 0,t'1) approaches a finite non-zero limit as \t'1\ 
approaches (large b), eq. (0) implies dN^/db oc b~^ for large b. 
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At any point along the line of sight where the first derivative vanishes, a large curvature 
(or absolute value of the second derivative) of random fields of cosmological interest, like 
the optical depth or density fluctuation, should correspond to a rare event (being correlated 
with high maxima, which are themselves rare) while a small value is common-place. The 

asymptote is then a non-trivial consequence of the expectation that small It^^'I's are not 
rare (in the sense that P[r^ = 0, t'I = 0] is a finite non-zero number). 

The above argument is not completely general, however, because the condition that 
P[r| = 0,t'1 = 0] is finite and non- vanishing is not the only possible quantitative mani- 
festation of the expectation that |t^'|'s are not rare. More generally, we can state the fol- 
lowing condition: that -P/[t^ = 0, /{t'I)] approaches a finite non-zero value as / approaches 
/{t'I = 0), where / is a smooth analytic function of r'l around t'I = and Pfij'j^ = 0, f)dT'j^df 
is now the probability that and / take the prescribed values. We allow the possibility that 
the old probability density -P(t^ = 0, t'I) vanishes at t'I = 0, but assume it is non-singular at 
that point. One can replace P in eq. (^ by Pfdf/dT'l and Taylor expand / around t'I = 
using / = Y^'^o fiT'l\ Eq. (|]) then implies that dN^/db oc b~"^ for sufficiently large b, with 
m > 5. For instance, suppose fi = for all i less than some number n (excluding /i=o, which 
is irrelevant for Pfdf/dT'l) and fi^O for i = n (with no additional condition imposed on all 
other /i's), then dNu/db oc fe"*" in the large b limit, with m = 5 + 2(ri — 1). Q 

The b~^ tail mentioned earlier corresponds to the case where /i is non-vanishing. We 
will show in §^ that this is the case for Gaussian random optical depth fluctuations. 

How large does b have to be for the asymptote to take over? One can define the 
following velocity scale: 



-'high 



\ ' dT'l 



p 



nfnT\ 



lln—l 
L 



(3) 



where all terms are evaluated at vanishing and t'I, and /j=„ is the first non- vanishing fi 
defined before. The quantity fehigh gives us an estimate of how large b has to be for the onset 
of the 6"™" asymptote. In other words, it tells us how low t'I has to be for P / {nfnT'U^~^) 
to be well-approximated by a simple non-vanishing constant. Note that in cases where the 
partial derivative with respect to t'I in eq. (|^) vanishes, for instance in the linear theory of 
analogous quantities can be defined using higher derivatives. 

How does the above theoretical prediction fare with observations? 



•^The analyticity condition on / is an important one which allows us to Taylor expand. Consider a 
counter-example where / = t£' with fc < 1, then df/dr^' diverges at t^' = 0, and our condition that 
P/[t^ = 0, fir'l — 0)] is finite would imply the old probability density P{t'^ — 0, t'I) diverges at r'l — 0. We 
take it to be a plausible assumption that P{t'^ = 0, r^') is non-singular. 
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From a sample of 790 Lya lines with 12.5 < log(A^Hi) < 14 taken from the published 
line-lists of Hu et al. ( |1995| ), we compute the cumulative b distribution and compare it with 
a single power-law dN/dh oc 6^*" using only lines above a value fe^ut using a Kolmogorov- 
Smirnov test ( [Press et al. 199'3 ), which produces the probability (pks) that a realization 
of the theoretical distribution would produce a dataset which is more disparate than the 
observed distribution is from the theoretical one. We note that the four lines of sights in the 
sample all have similar redshift ranges (see Table |l|). For lines with 6 > 30 km s~^(368 lines), 
the Pks > 0.005 range (at 0.01 increments) for m is [-4.12,-3.04]; for 6 > 40 km s~^(180 lines), 
it is [-5.03,-3.11]; and for 6 > 50 km s'^^M lines), it is [-6.10,-2.89]. We note that restricting 
the h range to higher values, one eventually can obtain a sample which is consistent with any 
power-law slope. Thus while the observed distribution is consistent with the — m (m > 5) 
power-law for sufficiently high 6's, it does not require one. Within the context of our theory, 
~ 40kms~^ can be regarded as a lower-limit to 6high defined in eq. (|^). 

Before we move on to the next section, an important caveat on the peak picture of an 
absorption line: in hierarchical clustering models, a given peak in the optical depth does not, 
in general, have an exact Voigt-profile shape. Standard profile fitting routines might fit the 
peak with a profile as described in eq. (0), together with a few smaller profiles to fill in the 
"wings" of the peak. If this occurs frequently enough, a prediction based on eq. (0) might fail 
to match the observed h distribution. However, the success of the peak picture in another 
context, the column density distribution, lends support to it. We will take it to be our 
working hypothesis, which can be checked through detailed comparisons with simulations. 
It should be noted also that higher order terms in the Taylor series expansion in eq. (H) 
contain information about departure from the Voigt-profile shape, and their statistics could 
in principle be computed. We leave it for future work. 



3. Thermally Broadened — or Not? 



The high velocity tail discussed in the last section implies that there always exist very 
broad absorption lines. An obvious question: should the high h value be taken to indicate 
high temperature, as is traditionally assumed? 

The answer is: not necessarily. In general, there inevitably exist absorption lines much 
broader than the thermal broadening width. Let us try to understand its physical origin. 
The optical depth r at a given velocity uq is given by (see |Hui et al. 1996|) : 



l + z 



du 
dx 



-1 



(Tadu , = aao- — ^exp[-(M - uqY /hx 



(4) 
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where nni is the proper number density of neutral hydrogen, z is the mean redshift of interest, 
X is the comoving spatial coordinate and the integration is done over the velocity u along the 
line of sight. The Jacobian \duldx\ multiplying the proper density rini gives us the neutral 
hydrogen density in velocity-space, and the summation is over multiple streams of x's at a 
given u. 

The thermal profile is given in the second equality, with (Jao being the Lyman-alpha cross 



section constant ( [Rybicki fc Lightman 19791) . The width of the profile is 6^ = yS/csT/mp 
where T is the temperature of the gas, fc^ is the Boltzmann constant and vny is the mass of 
a proton. Turbulence broadening can in principle be included simply by defining an effective 
temperature (T = Tturb. + ^thermal), and we will use the term thermal broadening to refer to 
both. 

It can be seen from eq. @ that if the HI number density in redshift space tt-hi \duldx\ is 
sharply peaked (e.g. a Dirac delta function), r has the shape of the thermal profile around the 
peak, and the measured h as proposed in eq. (|l|) would coincide with the width 6^. However, 
in the opposite limit in which nYa\dul dx\ is varying slowly around a peak, nYa\duldx\ can 
be taken out of the integral, and the result is that r does not have the shape of the thermal 
profile in general, and its width is determined by the scale of variation of nYi\\dul dx\^ which 
is larger than 6^ in this case. The measured b defined in eq. ([l|) then reflects the velocity 
structure of the peak in HI number density, rather than the temperature of the gas. The 
existence of absorption lines with high b (> 67-) is unavoidable as long as structure formation 
models allow fluctuations on large scales, in other words, low-curvature fluctuations. 

The arguments above set 6t as the lower limit to the observed width of an absorption 
line, but there always exist lines which are much broader. Moreover, if the HI number 
density in redshift space is intrinsically smooth on scales larger than 67-, most lines would in 
fact be wider than the thermal broadening width. Possible sources of such smoothing are 
the Jeans-smoothing and the effective numerical/observational resolution, to which we will 
return. 

In the next section, we develop a linear perturbative analysis of the optical depth fluctu- 
ation (small fluctuation limit), which exhibits all these possibilities, and allows us to predict 
the full b distribution, from narrow to broad lines. It will be demonstrated that the tem- 
perature of the intergalactic medium influences the overall b distribution not so much by 
giving a thermal width to every absorption line, but by determining the amount of smooth- 



*Press & Rybicki (1992) arrived at the conclusion that the thermal interpretation of the observed high b 
values is inconsistent with the hypothesis of heating and photoionization equilibrium by a UV background 
using a different argument: essentially baryon counting. 
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ing of the optical depth field, which affects the relative probability of high- curvature versus 
low-curvature fluctuations. We emphasize, however, that the key result in namely the 
existence of the b~"^ asymptote, is independent of the validity of the linear calculation. 



4. A Linear Analysis of the Optical Depth Fluctuation 

Assuming that t{u) = f[l + 6r{u)] where f is the mean optical depth and 6r ^ 1, the 
quantities and r'l in the expression for dNu/db in eq. (|^) can be replaced by 6'^ and 6'^ 
respectively. 

As we will show below, for cosmological models with Gaussian random initial conditions, 
the linear fluctuation 6r{u) is itself Gaussian random, and so the probability distribution 
P{S'^ = 0,6") is given by the product of two Gaussians, with 6'^ set to 0: {2n^ {5'^) {S"'^))^^ 
exp[— (5:5-'^/ {2{6"'^))], where () denotes ensemble averaging. The normalized b distribution then 
follows from eq. (0): 

b^^ 2 

bt^^^, (5) 



dN Abl 
-dh = ^^"P 



64 



where the normalization is chosen such that J^{dN/db)db = 1.Q 

The b~^ high b asymptote is clearly exhibited. As is shown in §^ this is a result of 
the fact that P{S'^ = 0,6"), which is proportional to exp[— 6^/6^] here, approaches a flnite 
constant in the large b limit. 

This linear Gaussian model also predicts the presence of a very sharp cut-off at low-6's; 
such a cutoff has been reported from observations (|Hu et al. 1995| ; |Lu et al. 1996| ; [Kirkman 



Tytler 1997| ; [Kim et al. 1997|) . This distribution, with fo^r = 26.3 km s ^, is shown 



super-imposed upon the observed distribution from one line of sight (QSO 0014-1-813) in 
Fig. |I[ Note how a single parameter b^ controls both where the low b cut-off and the high b 
asymptote take over. It is determined by the rms fluctuation amplitude of 6". Let us derive 
its dependence on cosmological parameters. 

Returning to eq. (ffl), let us flrst note that in the linear regime, one can ignore multi- 



''Thc non-normalized distribution, or the total number of absorption lines, is interesting in its own right. 
It turns out that to model it properly, at least one important selection effect has to be taken into account, 
namely typically only lines above a certain column density are included. The above formalism can be 
modified to accommodate it by including a threshold to Sr, which is correlated with 5". Such a modification 
would in fact alter the shape of the normalized b distribution as well. We leave complications due to this 
selection effect as well as others for a future paper. 
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pie streaming and therefore the summation. Ionization equihbrium imphes that the nm 
is determined by the local baryon overdensity, let us call it 6, and the local tempera- 
ture through rim oc [1 + (5]^T~°-'^. The temperature T is typically related to 5 by T = 
To(l + S)'^"'^, where Tq is the mean temperature at 5 = and 7 is determined by reioniza- 
tion history ( |Hui fc Gnedin 19961) . The velocity u is related to the spatial coordinate x by 
u = H{x — x)/{l + z) + fpec where x is the mean position of interest and H is the Hubble 
constant at redshift z, and fpec is the peculiar velocity along the line of sight. 

Collectively, the above relations imply that the optical depth fluctuation 6r is simply 
determined by two random fields: S and Vpec- The reader is referred to Hui et al. ( |1996| ) for 
details. Keeping only terms to first order in S and Vp^c, one obtains from eq. 



W{u-uo) 



[2 - 0.7(7 
1 



l)]S 



dv 



pec 



du 

exp[-(M-Mo)V&efr] 



+ (7 - 1) 



4 du^ 



W{u — uo)du 



(6) 



where 5r = {T — f)/T, and 6^0 = y'^ksTQ/mp is the thermal broadening width at temperature 
To, and W is simply a Gaussian smoothing window. The quantity fogfr should be set to 
strictly speaking, if one were to derive the above from eq. (^. However, we allow them to 
be different for the following reasons: 1. the observed spectrum is often convoluted with a 
Gaussian resolution window (strictly speaking the convolution is operated on exp(— r), but 
in linear theory the same convolution is applied to 6r) ; one can also model the resolution of 
the simulated spectrum similarly; 2. to relate the above quantity directly to cosmology, it is 
convenient to replace the baryon overdensity 6 by the dark matter overdensity (assuming the 
universe is dark matter dominated), which we will denote by S from now on, but a smoothing 
kernel has to be applied to the latter to take into account the smoothing of small scale baryon 
fluctuations due to finite gas pressure, which can also be approximated by a Gaussian ( [Hui| 
fc Gnedin 1996| ) (similarly for the Vpec field, which is related to 6 by fpec = —d'V~'^6/dx in 



linear theory, where the dot denotes differentiation with respect to conformal time). The 
combination of the three different Gaussian kernels, due to thermal broadening, resolution of 
observation/simulation and smoothing by finite gas pressure, is itself a Gaussian with width 
bes given by: 



■>To 



hi + 6L + h] 



13 kms 



lO^K 



where 



; ^res 



FWHM 



2Vln2 



6j = 24kms-Vj 



7 



1.5 



lO^K 



(7) 



where 671, is the thermal broadening width given before, fej-es is the resolution width which is 
related to the commonly quoted resolution FWHM by the factor given above, and 6j is the 
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baryon smoothing scale. If the factor /j were set to 1, 6j is exactly the Jeans scale in the 
appropriate velocity units (6j = 2fjHkJ^/{l + z) in the convention of [Gnedin fc Hui 1997| ). 
However, as shown by Gnedin & Hui ( |1997D , the correct linear smoothing scale is in general 
smaller than, but of the order of, the Jeans scale, if one waits for sufficiently long after 
reionization. For example, if reionization occurs at z = 7, /j ~ 0.5 at z = 3. It is interesting 
to note that both fe^Q and 6j scale as the square root of the temperature. 



, assume that both fields 8 and v^^c 



To establish the Gaussianity of 6r{u) in eq. 
are Gaussian random in the three-dimensional comoving space x, which is expected for a 
large class of inflationary cosmological models. One dimensional projection x to x preserves 
Gaussianity, and so does the comoving space to redshift-space mapping, to the lowest order 
(i.e. one can equate a change of variable from u to a; or vice versa with the simple linear 
transformation u = H{x — x)/{l + z), if one is only keeping flrst order terms as in eq. 

0) . These two facts mean that both 6 and fpec are Gaussian random in u space. Finally, 
all operations on 6{u) and Vpec{u) in eq. @, differentiation, addition and the Gaussian 
convolution, are linear which means 6r{u) itself is Gaussian random. 

Therefore, simply put, 6r is a Gaussian random fleld which is equal to [2 — 0.7(7 ~ 

1) ]6 — t;' + (7 — l)6f|j(5"/4, smoothed on the scale of bcs given by eq. (0). Similarly 6'^ is 
itself Gaussian random and its statistical properties are completely specifled by its two-point 
function. 

It is straightforward to compute {S'^"^), which gives in eq. (|^): 



where 



and 



a 



[(1+^)1 


4 


H 




- 0.7(7 - 


-1) 



-a 



2 1 



2 2 1 ■ 



d\n.a 



0.6 



U2 7 



l + zf 



/■oo 

/ A^0+^P(k) exp 

JO 



4 

y/2H 



H 



off 



6cff(l + z) 



(9) 



(10) 



A few symbols require explanation: is the linear growth factor, a is the Hubble scale 
factor, P{k) (to be distinguished from P[r|^, t'I] used earlier) is the linear-extrapolated power 
spectrum today as a function of the comoving wavenumber /c, VLm is the matter density, and 
redshift dependent quantities such as H and VL^n are to be evaluated at z. 

Given a cosmological model together with a reionization history, can be calculated 
from flrst principle, and eq. (|^) gives us the linear prediction for dN/dh. Let us see how 
the prediction fares with observations, and then sort out the somewhat complicated (multi-) 
parameter dependence exhibited above. 
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4.1. Comparison with the Observed Distribution 

We compare the hnear theory prediction for the differential b distribution in eq. 
with datasets from six QSO lines of sight from three different studies ( [Hu et al. 1995| ; |Ei: 



et al. 1996| ; [Kirkman fc Tytler 19971) . The theoretical prediction is compared individually 
with each line of sight, because in principle, the parameter can evolve with redshift. In 
Table ^ we give the range of b^- with pxs >1%, testing models with b„ in the range 10-30 
km s~^, at 0.1 km s~^ increments. Ranges were found for the four QSO datasets of Hu 
et al. (|1995|) , with acceptable b^ values in the range 23.3-28.8 km s~^. Of these four datasets, 
three produced a range of b„ with pxs > 10%. Data from the QSO HS 1946+7658 line-of- 
sight ( Kirkman fc Tytler 1997 ) and from the QOOOO-26 line-of-sight ( |Lu et al. 1996| ) were not 



compatible with with the theoretical distribution for the range of 6o-'s tested, with maximum 
probabilities of pxs = 4 x 10~^ and 3 x 10~^ at b^j = 26.1 and 23.7 km s~^, respectively. The 
differences found between the above studies could be partly a result of different profile-fitting 
algorithms. 

However, if we assume that b^ does not evolve significantly over the redshift range 
covered by the Hu et al. (|1995|) datasets (they are at similar, but not exactly the same, 
redshifts), grouping them together diminishes the agreement, and the maximum pxs reduces 
to 5 X 10~^. We interpret the above findings mainly as an indication that higher order 
corrections ignored in the linear analysis are not negligible, as we will demonstrate in §ff]2 
for a reasonable cosmological model. 

Nonetheless, that a single value parametrization of the b distribution can fit some of the 
datasets at all is interesting (as opposed to three parameters in the case of the commonly 
used truncated Gaussian). Let us also emphasize that the validity of the peak picture in 
general, and the high b asymptote in particular, as discussed in §0, are independent of the 
validity of the linear perturbative expansion. 



4.2. Physical Interpretation 

To quantify the departure from linearity, one can compute (5^) which is the rms 
fiuctuation amplitude of the optical depth 6r. For the CDM model with as = 0.7, Qm = 1, 
a present Hubble constant of 50kms~^ Mpc~^, and assuming b^s = 18kms~^ (for resolution 



FWHM of about 8kms-\ Tq = WK, 7 = 1.5 and /j = 0.5 in eq. 0), ^{5^) = 3.6 at 
redshift z = 3. The above parameters are known to give a column density distribution that 
agrees with observations (e.g. [Hernquist et al. 1996| , [Zhang et al. 1996| and Pui et al. 1996| ). 
It is interesting to note that for the same parameters, the rms fiuctuation amplitude of the 



- 11 - 



matter overdensity S is 1.8, a factor of 2 smaller: large scale coherent peculiar velocity flow 
helps enhance the optical depth fluctuation. 

Although the linear theory calculation fails, it may still help us understand the broad- 
stroke parameter dependence of the b distribution. The key parameter in the calculation is 
{S"'^), which gives an indication of how large the average b is expected to be, and how much 
spread one should expect in the b distribution (eq. 

Let us dissect the parameter dependence of {6""^) in eq. (||). Roughly speaking, af is of 
the order of (Jok'^ff with kes oc H/[bes{l + z)] (see eq. |TD|). Eq. (W) and (B) then imply 



assuming that a, fn and bro/bes are of the order of unity. For the same CDM model discussed 
above, an exact calculation yields b^ = 14kms~^, which is rather low compared with the 
observed b„ in the range 23-29 km s~^. We will return to this point below. 

The Smoothing Scale 

Clearly, the one important parameter determining the relative number of high b to low 
b absorption lines is the smoothing scale feeff, which is the only relevant velocity scale in the 
problem. The larger bcs is, the broader the lines are in general. Note that increasing bcs 
(decreasing fcefr) decreases ctq (eq. ll^l), which can only enhance the effect. 

The effective smoothing scale 6eff has contributions from three sources: thermal broaden- 
ing (including local turbulence) (&To)) observation/simulation resolution (6res) and smoothing 
due to finite gas pressure (~ Jeans scale; 6j) (eq. 0]). We can see from eq. (|^) the baryon 
smoothing scale 6j is in fact of the same order as for /j ~ 0.5, unless turbulence is 
important (in which case the effective Tq going into should be higher than the thermal 
To going into 6j). The precise value of /j depends very much on reionization history, the 
overall trend being the closer the epoch of reionization is to the time of interest, the smaller 
/j would be (see [Gnedin fc Hui 199^ ). At a redshift of 3, if reionization occurs before about 
redshift of 7, 6to is interestingly close to 6j. 

The resolution width bj-^s, on the other hand, can vary by large amounts. High quality 
Keck spectra have FWHM of 8kms~^, which corresponds to fores ~ 5kms~^ i.e. the resolu- 
tion effect is expected to be subdominant compared to the other two sources of smoothing. 
However, current simulations have effective resolution spanning about an order of magni- 
tude: fores ~ 3kms~^ — 30kms~^ ( |Cen et al. 1994| , [Hernquist et al. 1996| , [Zhang et al. 1996 



see, in particular, discussions in [Bond fc Wadsley 19971 ), which means for some of them, fo^ 
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could be the dominating influence on the average width of an absorption hne. This might 
explain some of the differences seen in their predicted b distributions ( pave et al. 1997| & 
[Zhang et al. 199^ ), although differences in their Voigt-proflle fltting techniques and reion- 
ization histories certainly contribute. One should also bear in mind that, in the presence of 
nonlinear corrections, the proper smoothing scale would no longer be exactly bes (even the 
Gaussian form of the smoothing would be changed), but the overall trend that higher b^^ 
and/or feres and/or bj imply broader lines is expected to remain true. 

In the context of the linear theory, one can see that if bes is dominated by br^ , ba would 
be of the order of (eq. [0], assuming -D+ctq ~ 1). This means the bulk of the absorption 
lines would have widths close to bx^, but the existence of much broader lines is also inevitable 
(eq. 0). Moreover, if 6efr is dominated by the baryon-smoothing scale or resolution, most 
of the lines would end up having widths larger than 6^0- 

The temperature of the intergalactic medium should not be understood as influencing 
the overall b distribution by giving a thermal width to every individual absorption line. 
Instead, it affects the b distribution by determining the amount of smoothing of the optical 
depth fleld (through 6^0 and which controls the relative probability of high- curvature 
versus low-curvature fluctuations (a higher Tq favors the later, thereby shifting the overall 
distribution to higher 6's; see eq. 0, and ||TT[|). 



Normalization of the Linear Power Spectrum 



The other (dimensionless) parameter influencing the b distribution in the linear theory 
is -D_|_(To (eq. ^^), which is simply the rms linear overdensity fluctuation at the redshift of 
interest. The more linear power a given model has at the smoothing scale 6efT, the lower 
bcr is, which implies on the whole narrower absorption lines. In other words, high-curvature 
fluctuations are more probable for a model with larger linear power. Whether this trend 
remains true after nonlinear corrections kick in is unclear without a detailed calculation. 
From the discussion in §^ we know that b = y^2/|(lnr)"|, independent of the validity of the 
linear expansion. One can take the magnitude of ((Inr")^) as an indication of how narrow 
the absorption lines are. Expanding it beyond the lowest order, subject to the extremum 
condition r' = 0, we have ((Inr")^) = (5^'^) — 2{6"'^6r)- The higher order corrections can in 
fact make ((Inr")^) smaller than the linear theory prediction. This is consistent with the 
fact that the linear prediction for the CDM model discussed earlier is b^ = 14kms~^ which 
is rather lower than the observed values (~ 23-29 km s~^; one could of course argue the 
whole perturbative expansion breaks down for the above model, but it might work better 
for models with less small-scale power such as the mixed dark matter models). The overall 
effect of raising the normalization of the linear power spectrum is therefore unclear. For the 
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same reason, it is difficult to make any qualitative statement regarding the redsliift evolution 
of the normalized b distribution, because it is unclear in what form the combination -D+ctq 
would appear in the expression for dN/ db if nonlinear corrections are included. On the other 
hand, one other source of redshift evolution is the evolution of the effective smoothing scale 
feeflf discussed earlier, which is mainly determined by the evolution of the equation of state of 
the intergalactic medium i.e. reionization history. 



Assuming fixed 6efr, the parameter Qrn enters into the linear prediction for the b distri- 
bution in two different places: and /o (eq. and p). Lowering has two effects: it 
decreases the growth rate D_^. and makes fn, which quantifies how the comoving-to-redshift- 
space mapping enhances fluctuations, smaller. However, for the same reason as discussed in 
the case of the power spectrum normalization, the qualitative effect of nonlinear corrections 
is difficult to guess without a detailed calculation. But the linear calculation suggests that 
the flrn dependence of the b distribution is inevitably tangled with the dependence on power 
spectrum normalization. One should also keep in mind that at sufficiently high redshifts, 
flm is close to 1 whatever is today. However, for open universe models, with flrn of 0.3 
today, Qm = 0.6 even at a redshift of 3. 

J 

Lastly, one parameter that is curiously missing from the linear theory calculation is the 
ionization background J. It is easy to see that, quite generally in fact, since one is only 
interested in the shape of r around an absorption peak, the overall normalization of J does 
not play a role in determining the b distribution. This is true, however, only in so far as the 
range of J being considered is small enough that it does not affect severely the selection of 
absorption lines into one's sample (e.g. because of finite signal to noise, etc). 

5. Discussion 

We develop in this paper a theory relating the 6-parameter of a Lja absorption line to 
the curvature of a peak in the optical depth fluctuation. It is shown in under the condition 
that low-curvature optical depth fluctuations are common, an asymptote of dN/ db oc b~"^ is 
in general expected in the broad-line limit, with m > 5. This is independent of details of 
the cosmological model, or the smallness of the fluctuation. The observed b distribution is 
consistent with the onset of the b~^ asymptote for b larger than fohigh ~ 40kms~^ (eq. 0]). 
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A perturbative analysis of the optical depth fluctuation is developed, that, in addition 
to illustrating the above asymptote, predicts a sharp low b cut-off. This is a reflection of 
the fact that high-curvature fluctuations are strongly suppressed in a Gaussian random fleld. 
Of the six datasets that we compared with the linear theoretical distribution, four were 
consistent with the theory (with pxs > 1%) for a range of b^. The other two were consistent, 
at best, at the pxs ~ 0.3% level. The fact that a single value parametrization as predicted by 
the linear analysis, as opposed to the commonly used three-parameter truncated Gaussian, 
works for some of them is interesting. However, we interpret the general disagreement as an 
indication of non-negligible nonlinear corrections. 

The effective smoothing scale feefr, which arises from a combination of thermal broaden- 
ing (including turbulence), baryon smoothing due to flnite gas pressure and flnite observa- 
tion/simulation resolution (6^0 5 ^J, ^res in eq. 0), is the key parameter that determines the 
overall width of the absorption lines. The higher fogfr is, the broader the absorption lines are 
in general. The baryon smoothing is shown to be at least as important as thermal broaden- 
ing for sufficiently early reionization, unless turbulence is a significant source of broadening. 
If reionization occurs very close to the redshift of interest, however, fe^o could dominate, 
because it takes a while for bj to grow to an appreciable value after reionization ( pnedin 
lui 199 



The importance of fej-es varies a lot between different observations/simulations. The 
resolution effect is probably unimportant for high quality Keck data, while some of the 
current hydrodynamic simulations might produce resolution-dominated b distributions. The 
fact that Keck-quality data are not resolution-dominated also means one can potentially use 
the observed b distribution to probe reionization history through its dependence on 6^0 and 
bj ( |Haehnelt fc Steinmetz 1997| , see also [Hui fc Gnedin 1996| ). 



As we argue in while the thermal broadening width does provide a lower limit to the 
width of an absorption line, the existence of much broader lines is an inevitable consequence 
of the fact that low-curvature fluctuations are statistically favored. The simple thermal 
interpretation of all the observed b values is therefore not viable, at least in the context of 
current structure formation models. The temperature of the intergalactic medium affects the 
overall b distribution not by giving a thermal width to every individual absorption line, but 
by determining the amount of smoothing of the optical depth fluctuation. Any attempt to 
understand the full b distribution of the Lya forest and its dependence on reionization history 
must involve understanding the statistics of the optical depth fluctuation, and its relation 
to the underlying cosmological model. A flrst attempt has been made here to explore how 
cosmological parameters such as flm and the power spectrum influence the b distribution 
in the context of linear theory. Further research taking into account nonlinear effects is 
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worth pursuing. The velocity structure of absorption peaks, quantified here by the b value, 
may provide important constraints on reionization history, as well as other cosmological 
parameters. 
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Fig. 1. — Comparison of the linear theory prediction for the differential 6-distribution, 
using = 26.3 km sec~^ (eq. p ) with the observed distribution of 194 Lja lines with 
1012.5 < jY^j < lO^^cm-2, from Q0014+813 (Hu et al.l995; see Table U) . The theoretical 
distribution qualitatively and quantitatively matches the data, with a cut-off at low b values, 
a quasi- Gaussian peak, and a high-6 tail. Statistically, the fit is acceptable (probes = 0.16). 



This preprint was prepared with the AAS M5t;X macros v4.0. 
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Table 1. b Data-sets 



QSO 


(•^lines) 


No. of b values 


bcr (km s ^) 


Max(pKs) 


6,(Max(pKs)) 








[Pks > 1% range] 




(km s--*^) 


HS 1946+7658(A) 


2.7 


328 




4 X 10-3 


26.1 


Q0014+813(B) 


2.95 


194 


[25.5, 27.2] 


0.17 


26.3 


Q0302-003(B) 


2.86 


198 


[25.7, 28.8] 


0.25 


27.5 


Q0636+680(B) 


2.75 


219 


[23.3, 24.8] 


0.09 


24.2 


Q0956+122(B) 


2.85 


176 


[24.7, 27.2] 


0.20 


26.3 


Q0000-26(C) 


3.7 


287 




3 X 10-3 


23.7 



References. — (A) [Kirkman fc Tytler 1997| : (B) |Hu et al. 1995| : (C) |Lu et al. 1996 
Only lines with 10^^'^ < A^hi < 10^^ cm-^ are included in our samples. 



